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1 Introduction 



Short-distance phenomena in strong-interaction physics are usually described 
in the framework of a perturbation expansion. The need for high precision is 
one reason why we may want to reach beyond the limitations of this framework. 
But there are also interesting questions regarding the structure and meaning 
of the perturbation series itself which naturally lead us to investigate long- 
distance, non-perturbative aspects of hard processes in QCD. The purpose 
of this article is to exhibit this relation of perturbative and non-perturbative 
physics, and the kind of phenomenology it has entailed. Some of the details 
which we have to skip here can be found elsewhereu 

To state this more clearly we consider an example from deep-inelastic 
neutrino-nucleon scattering. In the parton model the Gross-Llewellyn-Smith 
(GLS) sum rule for the structure function F3 expresses the fact that the proton 
consists of three valence quarks. In QCD there are finite corrections to this 
sum rule, so that 



I GJ * = \j\ (F» p +F^)(x,Q) 



«• _ 3.58 f^) 2 - 18.98 - . . . - % + 0(1/Q 4 

TT \ TT / \ TT / Q Z 



,(1) 



where a s — a s (Q) and the leading higher-twist correction C2 has been esti- 
mated to be 0.1 GeV 2 . Conventional parlance would say that the sum rule 
has a perturbative contribution (the series in o: s /tt) and a non-perturbative 
one, but this is an imprecise characterization of Eq. (^), because the pertur- 
bative series diverges for any value of a B . What is the numerical value of the 
"perturbative contribution" ? 

There are several reasons for why the series expansion might diverge. 
The divergence^ that has most implications is known as the infrared (IR) 
renormalonfltrel (The other known sources of divergence are: instanton-anti- 
instanton pairs, which are suppressed in QCD; ultraviolet renormalons, which 
- at least in principle - can be disposed of for practical purposes. We briefly 
discuss ultraviolet renormalons in Sect. ||) For the GLS sum rule it can be 
shown that the coefficients r n of diverge (due to IR renormalons) for 

large n as 

r n o,{^y n \n\ (2) 

where p — 2, /3q = (11 — 2iV//3)/(47r), Nf the number of massless quarks, and 
b some constant which we assume here to be zero for simplicity. A divergent 
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series expansion is a useful approximation if it is asymptotic to the quantity 
which it represents. QCD perturbative expansions have never been proven to 
be asymptotic, but it is a good idea to proceed with this assumption on good 
faith. Assuming furthermore that 



with K n oc r n the best approximation occurs at N = No ~ \p\/(2p$a s ) and 



With p = 2 this is of the same order of magnitude as the power correction 
C2 / Q 2 in Eq. ([[]) . Such power corrections are referred to as "non-perturbative" , 
because they are exponentially small in the strong coupling a s (Q). Without 
a summation prescription, however, the numerical value of the "perturbative 
contribution" is not unique, although it is determined to an accuracy (A/Q) 2 , 
which is small, when Q is large. It therefore seems that to include these power 
corrections consistently we only need to figure out the correct prescription 
to sum the perturbative series and add to the sum the higher-twist correc- 
tion C2/Q 2 - but this is wrong! Defining what C2 is requires great care, and 
once the power correction is properly defined, the summation prescription for 
the perturbative series is implied by this definition.LrEm Two important con- 
sequences follow: the perturbative series and power ( "non-perturbative" ) cor- 
rections are not independently defined, they are related; if the series diverges 
as in Eq. (ph, then there must be a power correction of order (A/Q) p , whose 
precise definition fixes the ambiguity in defining the perturbative expansion. 
We can use this to obtain some insight into power corrections using nothing 
but the rules of perturbative QCD. 

This becomes much clearer, if we take into account the physics origin of IR 
renormalons. The sum of Feynman amplitudes at a given order in perturbation 
theory, which give the perturbative expansion of the GLS sum rule, is IR 
finite and depends only on the large scale Q. On dimensional grounds the 
average loop momenta must therefore scale with Q. However, the coefficient 
of proportionality may depend strongly on the order of perturbation theory. 
Suppose in ann+1 loop contribution to /qls we have integrated over all loop 
momenta but one momentum k and that the result of the n loop integrations 
is proportional to (3q \n' l (Q 2 /k 2 ). The coefficient (3q appears unmotivated at 
this stage and will be explained in Sect. ||, but for now we may only note that 
the scale dependence of the coupling is an obvious source of logarithms. Then 




(3) 
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the dominant contributions to the final integral over k come from k 3> Q and 
k -C Q, because of the large logarithmic enhancements in these regions. (The 
contribution from large k is related to ultraviolet renormalons and we ignore 
it in the following.) If I(k), the integrand for the final loop integration, goes 
as I(k) ~ k p ~ 4 for k <C Q, we find 

d 4 fc/(fc)/3£ln"(Q7fc 2 ) ~ ) n! (5) 



P 

as in Eq. (g) with typical k ~ Qe~ n / p . The contribution to r„a" +1 from 
fc < few- A is again of order (A/Q) p . The rules of perturbative QCD cannot be 
assumed to account for this small loop momentum region. This is not a prob- 
lem in conventional 1-loop or 2-loop calculations, since the power-suppressed 
contribution from small momenta is much smaller than the dominant con- 
tribution of order a s or a 2 s . But since we are now interested in such small 
power-suppressed effects, we should better consider the small loop momen- 
tum contributions as part of the low-energy matrix elements of higher-twist 
operators in the operator product expansion (OPE) of the GLS sum rule. 

We therefore define Ci as the matrix element of the relevant twist-4 opera- 
tor (strictly speaking, as the product of the matrix element and a coefficient 
function) which includes all low momentum contributions with k < /i for some 
/j <C Q, but larger than A. The contribution from k < fi has to be accordingly 
subtracted from the perturbative expansion. Since the higher-twist operator is 
quadratically ultraviolet divergent, we know that C2 = c/i 2 + 0(A 2 ) for /x ^> A 
and hence we can rewrite Eq. (Q) as 



GLS 



i-E^" +1 + ^}-§ + o(i/q 4 ) 



(6) 



with the subtraction of low momentum regions taken into account in the curly 
brackets. For \i sufficiently large compared to A, c itself has a perturbative 
expansion in a s . This expansion is exactly such that in large order it can- 
cels the IR renormalon divergence of the coefficients r„ so that the two ex- 
pansions in curly brackets combined are convergent - more precisely, the re- 
maining IR renormalon divergence causes a summation ambiguity of higher 
order in the 1/Q expansion. By defining C2 accurately, we succeeded in sum- 
ming the divergent series by eliminating the divergence altogether! However, 
the short-distance (k > //, "perturbative") and long-distance (k < ^1, "non- 
perturbative" ) terms in Eq. ([|) depend separately on fi and remain related 
just as before. But we now see that it is really the small-fc behavior of the 
1-loop integrand I{k) rather than factorial divergence which determines the 
magnitude of the power correction. 
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The analogy with conventional scale dependence may help understanding 
why the IR cutoff dependence of Feynman amplitudes can give information 
on non-perturbative power corrections. The former scale dependence relates 
different orders in perturbation theory and is often used to estimate the size of 
unknown higher-order corrections; the /i-dependence in Eq. (^) relates differ- 
ent orders in the 1/Q expansion and may be used to estimate the size of the 
first power corrections. These estimates are parametrically correct, but quan- 
titatively only indicative. The important point is that if the scale dependence 
is large, then so must be the higher-order correction in a s , or 1/Q. (It is often 
stated that the perturbative expansions are scale-independent to all orders of 
perturbation theory. Because the expansion is divergent this is only formally 
true. Any attempt to interpret the series numerically introduces the kind of 
scale dependence or prescription dependence exhibited in Eq. (||).) 

There exist several ways of making use of this connection between IR 
renormalons and power corrections: 

- formal: if power corrections can be analyzed with operator product ex- 
pansion methods, the saajf methods can be used to determine the IR 
renormalon divergence. BoO This goes as far as fixing all parameters of 
Eq. (^|), including subleading corrections in 1/n, but excepting the over- 
all constant of proportionality, since most of the structure of Eq. (^) is 
determined by logarithms which can be controlled by OPE and renor- 
malization group methods. This is perhaps of less interest phenomeno- 
logically, except to remind us that combining perturbative expansions 
with higher-order terms in the OPE is subtle unless we can argue that 
the matrix elements of higher-dimension operators are much larger than 
the low momentum contributions in perturbative Feynman amplitudes. 

- qualitative/scaling: here we begin with the divergence of the perturbative 
expansion and deduce from it the scaling with Q of power corrections. 
Some power corrections can be missed in this way, but usually there is an 
identifiable reason for this. The great advantage of this method is that 
the quantity does not have to admit an operator product expansion, the 
only requirement being that it has a short-distance scale, i.e. a perturba- 
tive expansion to begin with. In general this is a poor substitute for a full 
understanding of power corrections in terms of operators, but in some 
cases this is the only method known to this date. Rather than speaking 
of divergent perturbative expansions, we could directly investigate the 
small momentum behavior of Feynman amplitudes. This makes appar- 
ent the close relation of this approach with the methods of perturbative 
infrared factorization, but now extending this notion beyond the study 
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of logarithmic collinear and soft infrared sensitivity at the leading power 
in the hard scale. The Feynman amplitudes are implicitly presumed to 
indicate the correct scaling behavior of non-perturbative corrections. 

- quantitative: this is the most interesting, but also most delicate of all ap- 
plications. From what has been said it seems impossible to obtain quan- 
titative information on power corrections by perturbative methods. We 
cannot compute the small momentum contribution of arbitrarily compli- 
cated Feynman amplitudes, and even if we could, this would be of no use 
since this does not give the correct non-perturbative result. However, we 
can imagine a situation in which the subtraction term in curly brackets 
in Eq. (g) cancels almost completely all higher-order terms in the origi- 
nal pertubative expansion (say, for n > 1) for some value of /x, because 
the higher-order terms are already dominated by small loop momentum. 
If the value of n at which this occurs (say, fi = 1 GeV) is numerically 
large compared to A, the power correction C 2 /Q 2 in Eq. (g) might be nu- 
merically dominated by the first term in the expansion of cfi 2 in a s (/i) . 
In this case Eq. (g) is well approximated by a perturbative expansion 
truncated at n = 1 together with a power correction whose numerical 
coefficient is predicted by IR renormalons or perturbative IR contribu- 
tions! This is of course a rather idealized situation but we shall see that 
this logic provides an explanation of the sometimes puzzling success of 
models of power corrections based on perturbative infrared sensitivity. 
Since the power correction is really a parametrization of perturbative 
contributions, though originating at scales much smaller than the hard 
scale Q, we should expect that the "non-perturbative" power correction 
decreases as more terms are added to the perturbative expansion. 

The outline of this review is as follows: 

In Sect. |^ we show how the IR renormalon divergence can be character- 
ized with OPE methods. This is then illustrated by computing a set of fermion 
bubble diagrams, though in general we shall try to free the notion of renor- 
malons from its association with this rather special set of diagrams. The step 
to quantities without an OPE is made in this section, appealing to the more 
general concept of perturbative infrared sensitivity as already discussed in this 
introduction. 

Section |^ concentrates on phenomenological applications. In our opinion 
ideas based on, or inspired by, IR renormalons have had the most important 
impact on our understanding of power corrections to hadronic event shape 
measures; on perturbative effects in heavy quark decays and production due 
to the clarification of the role of heavy quark mass definitions; on the modelling 
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of twist-4 corrections in deep-inelastic nucleon structure functions. These three 
applications are discussed in some detail. Others can only be briefly summa- 
rized. 

Section ^ treats rather briefly the physics of ultraviolet renormalons. We 
conclude in Sect. ||. 



2 Infrared Renormalons - Basic Concepts 

2.1 The Borel Plane 

Renormalon divergence is often discussed in terms of the corresponding singu- 
larities in the Borel plane. We briefly introduce the relevant concepts. 

Given a quantity R and its series expansion, we define the Borel transform 
B[R](t) of the series by 



j-n 

1 => B[R]{t)=Y j r n - v (7) 



n=0 n=0 



If B[R)(t) has no singularities for real positive t and does not increase too 
rapidly at positive infinity, we can define the Borel integral (a s positive) as 

/>oo 

R= / dte- t/a ° B[R](t), (8) 



Jo 

which has the same series expansion as R. In QCD the Borel integral does not 
exist, since IR renormalons generate singularities on the integration contour. 
The non-existence of the integral is not a serious concern, however, since we 
would not have expected the Borel integral to equal R anyway, because of non- 
perturbative, power-suppressed effects. From Eq. (|§|) we see that the effect of 
singularities at finite t on R is also power-suppressed. Unless & is a negative 
integer, the correspondence between factorial divergence and singularities is as 
follows: 

r„ = KaT(n + 1 + b) <=> B[R](t) = (9) 

Because of this the divergent behavior of the original series is encoded in the 
singularities of its Borel transform. Hence, divergent behavior is often referred 
to through poles/singularities in the Borel plane. This language is particularly 
convenient for subleading divergent behavior. Note that larger a, i.e. faster 
divergence, leads to singularities closer to the origin t = of the Borel plane. 
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UV renormalons 

at t = — m//3o, m = 1,2, 



Instanton-anti-instanton 
singularities at t = 47r, 8ir, 



IR renormalons 

at t = m//3o, m = 2, 3, 



Figure 1: Singularities in the Borel plane of I1(Q 2 ). The singular points are shown, 
but not the cuts attached to each of them. 



To illustrate these concepts and survey the known singularities, we consider 
the correlation functions of two vector currents = qj^q of massless quarks 



(-<) / d i xe- i " x (0\T(j^x)j v (0))\0) = (q^-q 2 9 ^) U(Q 2 ) 



(10) 



with Q 2 = —q 2 . The singularities in the Borel plane are shown in Fig. |ll 
How general is this picture? In QCD, ultraviolet renormalon and instanton- 
antiinstanton singularities always occur at the locations indicated in the figure, 
independent of the particular observable, because the physics they represent 
(ultraviolet behavior and vacuum structure, respectively) is universally the 
same. The location of IR renormalon singularities depends on the specific 
observable. For observables derived from off-shell correlation functions, such 
as H(Q 2 ) and the GLS sum rule, IR renormalon singularities occur at integer 
multiples of 1//3q> because the OPE is an expansion in even powers of the 
hard scale. For observables derived from on-shell correlation functions, one 
can have power corrections suppressed by only one power of the hard scale. 
These would be related to a singularity at t = 1/(2/3q). In general, one csua 
construct infrared finite observables, which are arbitrarily infrared-sensitiveJiSI 
IR renormalons can then occur at any t and arbitrarily close to 0. 



2.2 IR Renormalons and the Operator Product Expansion 

From Eq. (||) we learned that factorial divergence arises if the typical loop 
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momentum k — ► as the order of perturbation theory increases. We take 
this to be the definition of what we mean by "IR renormalon" , i.e. an IR 
renormalon is a singularity of the Borel transform that is eliminated if all loop 
momenta are restricted to k > A no matter how small A is. By definition the 
analysis of IR renormalons is the analysis of small momentum contributions 
to Feynman amplitudes. This appears complicated since IR renormalons refer 
to large orders in perturbation theory and so it seems that we would need 
to investigate the infinite set of Feynman diagrams. We shall now see that 
the factorization properties of Green functions strongly constrain the form IR 
renormalon divergence can take.BQ 

For the remainder of this subsection we restrict ourselves to observables 
derived from off-shell (Euclidean) correlation functions; the function H(Q) de- 
fined in Eq. ( |l0| ) will serve as an example. The factorization of sma^lojao 
momentum regions is done for us by the operator product expansiorJlj'PH 13 ! 
(OPE), which for U{Q) reads 

n(Q) = C (a s , Q/n) + -L C d (a s , Q/fi) (— G 2 )( M ) + 0(1/Q 6 ), (11) 

where d = 4 and the leading power correction is given by the gluon condensate. 
The conventional way of interpreting this formula is to take Co as the series 
computed with the standard rules of QCD perturbation theory. This includes 
the small loop momentum regions which give rise to IR renormalons. It is 
conceptually more satisfactory to include these regions into the definition of 
the vacuum condensates - the OPE guarantees that this can always be done. 
The series expansion of Co then differs from the standard one and is free 
from IR renormalons. Both Co and the vacuum condensates are in this case 
separately well-defined. In the following we shall adhere to the conventional 
interpretation, but we will return to the second more appealing one later in 
this section. 

Since any loop momentum region with k < A can be absorbed into a series 
of vacuum expectation values, the IR renormalon contribution to Co must take 
the factorized form 

C R (a s , Q/ M ) = ij C d (a s , Q/fi) ^ d M(a s ) + 0(1/Q 6 ), (12) 

with the same coefficient function C d as in Eq. ( pd| ) and M(a s ) a (dimen- 
sionless) perturbative series independent of Q. The coefficient function of the 
gluon condensate, C d , satisfies the renormalization group equation 

{f^ +(3(a s )^--^^jC d (a s ,Q/^ = Q, (13) 
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where a s — a s (/i) and 7(013) is the anomalous dimension of the gluon conden- 
sate. Since Co and hence Cq R is independent of this implies that 

P(a s )^ + d -±pl)M {as ) = 0. (14) 

Here it is used that M being independent of Q and dimensionless cannot 
depend explicitly on jj,. The IR renormalon divergence must be contained 
entirely in M. We make the ansatz 

M(a s ) = u 7 : +1 K(a(3 ) n n! n h (l + — + <3(l/n 2 )) (15) 

n 

and insert it into Eq. (|l4|). We use (3(a s ) = —Po&l — ftictg — . . . and j(a s ) = 
7oa s + • ■ ■ and after shifting the summation index n in some terms and ex- 
panding in l/n, we obtain 

= £< +1 #(a/? ) r 

n 

I/O ' 1 ..,1 / " \ \ .....■) , .. 

(16) 





d 1 


nln b 






2 a 


Pi 


7o 


a 2 Pl 


2aP 



_--J Sl J + 0(lK) 

This implies that either K = 0, in which case there is no factorial divergence, 
or otherwise we must have 



2 dfh _ jo 

d' 2/3 2 2/3 



6 =77732 - — ■ ( 17 ) 



In the case of the gluon condensate we have d — 4 and 70 = 0, in which case 
we reproduce the result originally derived by Muelleru in a different way. Note 
that s\ and higher-order terms are also determined by the /3-function and the 
anomalous dimension .0 

It is important to stress the general nature of this result, since it ap- 
plies to all observables for which we know the structure of power corrections. 
There may be power corrections and no corresponding IR renormalons, i.e. 
K = 0, even if the higher-dimension operator has a power divergence (and 
hence a perturbative contribution) . This can occur in conformal theories such 
as supersymmetric Yang-Mills theory with four unbroken supersymmetries. 
However, if the divergence occurs, the n-dependence is completely determined 
by renormalization constants, theu-oecturbative expansion of Cd and the truly 



"non-perturbative" parameter K Jl3'&3 There is a common misconception that 
the involvement of the non-abelian /3-function coefficient /3q in the location of 
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the renormalon singularity is only a "conjecture" . This misconception is based 
on the set of fermion bubble diagrams (see the next subsection), which gives 
only part of 0q. It is indeed difficult to identify diagrammatically the factor (3$ , 
but the renormalization group allows us to bypass the difficulty. This should 
not be too surprising since the origin of renormalons is large logarithms. Hence 
we can prove the factor (3q , if the renormalon singularity exists, but we cannot 
rigorously prove that renormalons exist in QCD, for we cannot exclude that 
K = for some mysterious (hence entirely improbable) reason. 

Suppose we had excluded all loop momenta k < A from the coefficient 
functions. The IR renormalon would then reappear as the coefficient of an 
ultraviolet power divergence of the gluon condensate, although the condensate 
being defined non-perturbatively with the cutoff k < A, there is no need to sep- 
arate this divergence from the remaining condensate. It is however instructive 
to see how this works in an example. We cannot compute condensates analyt- 
ically in QCD, but the non^u^jJ^AD a-model in two space-time dimensions 
provides a nice toy model, t3CHlZl'Q't§ll3 which is solvable in a 1 /^-expansion. 
As QCD it has only massless particles in perturbation theory, but exhibits dy- 
namical mass generation non-perturbatively and a mass gap in the spectrum. 
It is asymptotically free, as is QCD, and m, the dynamical mass of the <r- 
particle, is the analogue of the QCD scale A. We cannot go into the details of 
the model here except to say that the vacuum expectation value of the square 
of the auxiliary field a{x), 



can be considered as the tr-model analogue of (a s G 2 ). Note that the restriction 
p 2 < fi 2 defim^the otherwise singular operator product a 2 . The integral can 
be evaluatedliZra with the result 

(a 2 )(^,m) = m 4 [Ei(ln A) + Ei(— In A) - In In A - ln(- In A) - 2 lE ] , (19) 

where je — 0.5772 ... is Euler's constant, Ei(— x) = — dter 1 Jt the expo- 
nential integral function and 



Note that F{x) = Ei(— x) — hxx has an essential singularity at x = but no 
discontinuity. 




l 



(18) 




(20) 
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Equation ( |19| ) is proportional to the dynamically generated scale to the 
fourth power, as expected, but for /i 3> m it develops power-like cutoff depen- 
dence. To see the emergence of renormalons, we expand the vacuum expec- 
tation value in powers of m/fJ, and the a- model coupling g(fj-) = —Pog(fi) — 
1/ ln(/i 2 /m 2 ). To perform the expansion we need the asymptotic expansion of 
F(x) at large x. For positive argument the asymptotic expansion is 



If the divergent series is understood as its Borel sum, the right hand side equals 
F. For negative, real argument, one obtains the asymptotic expansion 



Note the "ambiguous" imaginary part in the exponentially small term. The 
interpretation of Eq. (^2|) is as follows: the upper (lower) sign is to be taken, if 
the (non-Borel-summable!) divergent series is interpreted as the Borel integral 
in the upper (lower) complex plane. With this interpretation Eq. (^||) is exact 
and unambiguous. Inserting these expansions, the condensate is given by 



The expansion for large /i has quartic and quadratic terms in /i, parametri- 
cally larger than the "natural magnitude" of the condensate of order m 4 . The 
power terms inu arise from the quartic and quadratic divergence of the Feyn- 
man integral (|lq) , i.e. from loop momentum p ~ /i. The /i-dependence cancels 
with the /i-dcpcndcncc of the coefficient functions in the OPE. In particular 
the /i -term cancels with the coefficient function of the unit operator. The 
important point to note is that the condensate is unambiguous, but separat- 
ing the "perturbative part" of order /z 4 is not, since the asymptotic expansion 
for n/m ^> 1 leads to divergent, non-sign-alternating series expansions, which 
require a summation prescription. The "non-perturbative part" of order m 
depends on this prescription (via ±in in Eq. (|23|)). In a purely perturbative 
calculation, one would only obtain the divergent series expansion. The in- 
frared renormalon ambiguity of this expansion would lead us to correctly infer 




(21) 




(22) 
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Figure 2: The set of "bubble" diagrams consists of all diagrams with any number of 
fermion loops inserted into a single gluon line. 



the existence of a non-perturbative power correction of order m . For quan- 
tities without an OPE this is one of the main motivations for considering IR 
renormalon divergence. 

The renormalon ambiguity does not allow us in general to say much about 
the magnitude of the power correction which is determined by other terms, 
such as \n(2/g) in Eq. ([23]). The er-model is somewhat special in this respect, 
since the power-like ambiguities in defining perturbative expansions are also 
parametrically smaller in 1/N than the actual condensates. This tells us that 
some caution is necessary in identifying the magnitude of the "renormalon 
ambiguity" with the magnitude of power corrections. It is probably more 
appropriate to say that power corrections are expected to be at least as large as 
perturbative ambiguities. On the other hand, a similar parametric suppression 
of perturbative ambiguities does not seem to take place in QCD. 



2.3 The Large- flo Limit 

The best representative of renormalon divergence is the set oLfiermion "bubble 
diagrams" . Renormalons have originally been discoveredErau in this set of 
diagrams. It is important to bear in mind that the concept of renormalons is 
more general and that all diagrams eventually contribute to the overall constant 
K that appears for example in Eq. (|l5|). However, the bubble graphs are useful 
for explicit calculations. With a certain amount of extrapolation they have also 
turned out to give useful approximations to perturbative expansions in QCD. 

We consider again the current correlation function defined in Eq. (|To|), 
more precisely the Adler function D(Q) — 4ir 2 dH(Q)/dQ 2 , and compute the 
"bubble diagrams" shown in Fig. |2[ Any number of quark loops may be in- 
serted into the gluon line; each loop gives a factor (3 0f a s [ln(-fc 2 /^ 2 ) - 5/3], 
if the strong coupling is renormalized in the MS scheme, and (3of — —Nf/(6w) 
is the quark contribution to the 1-loop (3- function. i— in 

Bubble diagrams can be computed in several ways. Originallji23 : EJ the 
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Borcl transform was computed directly. The result is 



B[D]{u) = § (? e_5/3 



{-ifk 



l-(l-u) 2 ^ (k 2 ~ (1 -uff 



(24) 



with it = Poft. This has the singularities shown in Fig. [l] except for the presence 
of Pof rather than flo. The origin of the singularities is more evident, if we 
integrate over the loop momentum of the "large" quark loop in Fig. ^ and the 
angles of the gluon momentum k. Defining I — —k 2 /Q 2 , we obtain 



v r°° A] 

D =Y,o° 7 no 



Q 2 e' 5 / 3 



whereE3 

F(l) 

F(l) 



81 
3tt 

8/ 
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-f3 0f a s In I 



-]nl]l + (l + l) U 2 (-l) + ]nlbx(l + l) 



(25) 



/ < 1, 



3 4 1 + ln/ 



3 1 , A 1 

4 + 2 ln 0l 



+ (1 + LiaH-^-lniMl + r 1 ) 



I > 1 



(26) 



with 1^2(2;) the dilogarithm. The dominant contributions to the integral come 
from / 1 and I <C 1, because of the large logarithmic enhancements in 
these regions. There is a one-to-one correspondence between each term in the 
expansion for small (large) I and the IR (UV) renormalon poles in the Borel 
transform. For instance, the leading term at small I, 
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which corresponds to the IR renormalon pole at u = 2. 

We can verify explicitly the general result that this expansion can be in- 
terpreted as part of the gluon condensate. When k — > the gluon line in Fig. || 
can be "cut" , i.e. supposed to end in a slowly varying external field. The result 
of the computation is 27r/(3Q 4 ), the leading-order coefficient function of the 
gluon condensate. To verify that indeed 
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we compute .the (perturbative) gluon condensate in the bubble approximation 
and obtainEao 



(^G 2 )(fc< M ) = ^ Q J dee 



-f3ofa s In 
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= i^ 10/3 E«" +1 (^)"" ! - ( 3 °) 

n=0 ^ ' 

Combining this result with the coefficient function, there is agreement with 

The problem with the bubble approximation in QCD is that the renor- 
malon singularities are determined by the quark contribution to the /3-function. 
The general arguments tell us that (3of must get converted into (3q. We could 
add gluon and ghost bubbles as well, but this would still not give a complete 
result. In one way or another, recovering (3q in QCD leads beyond the approx- 
imation of a single dressed gluon line. Even if we could construct an effective 
charge analogous to QED, giving the complete (3q for every dressed gluon line, 
it is not clear what would be gained from such a construction, since the over- 
all normalization of renormalon divergence remains as elusive as before, and 
this constant is all that is not already determined by renormalization group 
arguments. 

H> esp>l"t e these diagrammatic difficulties it has been suggestedEj'C 2 ^! 2 ^ to 
use the quark bubble calculation with /3q / replaced by /3q as a realistic approx- 
imation to the full QCD perturbative expansion. Formally, this amounts to 
rewriting a perturbative coefficient at order as 

r n = r n0 + r nl N f + ... + r nn Nf = r [d n ffi + S n ] , (31) 

where d n — (~6n) n r nn /r , f3 = (11 — 2Nf/3)/(4n) and Nf is the number of 
massless quarks. The coefficients d n are then obtained from a calculation of 
fcrmion bubble graphs, while the remainder 5 n is neglected. For this reason, 
this procedure is often referred to as the "large-/?o approximation". It can 
also be viewed as an extension of Brodsky-Lepage-Mackenzie scale setting.tfl 
It is more an empirical observation than a well-reasoned statement that the 
remainder is indeed often found to be small compared to the f3o term (in the 
MS scheme). Since the large- /3q limit also incorporates the expected divergence 
of the QCD expansion, it has turned out to be a useful quantitative tool to 
estimate higher-order coefficients, in particular when the onset of divergence 
is rapid and when the observable depends only on a single scale. 



1G 



2.4 IR Renormalons and Power Corrections in Quantities without an Oper- 
ator Product Expansion 

Up to now we have been characterizing infrared renormalon divergence using 
known operator product expansion (OPE) methods. Much of the interest in 
IR renormalons during the recent 5 years comes from adopting a somewhat 
different point of view. Rather than concentrating on quantities whose OPE 
is well understood, let us take any "hard" quantity, i.e. a quantity that admits 
a perturbative calculation. If we are able to identify a particular pattern 
of IR renormalon divergence, we must replace the IR sensitive perturbative 
contributions by a non-perturbative parameter of similar magnitude. In this 
way, the suppression of power corrections can be obtained in a very general 
way. 

To our knowledge this possibility was first mentioned by Mueller ,0 who 
suggested to investigate hadronic event shape observables and jet cross sections 
from this perspective. The idea fell into oblivion at the time and it took another 
decade before renormalons in observables withoul-pj OPE were beginning to 
be studied in detalLiLrst in heavy quark physicsE2lEj and soon after for event 
shape observables. EH 

It is characteristic of off-shell processes that IR renormalons occur only at 
positive integer multiples of 1/(3q, which implies power corrections in powers 
of l/Q 2 and not powers of l/Q, where Q is the "hard" scale of the process. 
On-shell quantities have a large variety of infrared-sensitive regions in their 
Feynman amplitudes and the generic situation leads to IR renormalons at 
positive half-integers and integers and a series of power corrections in l/Q. 
Due to the variety of possibilities, we restrict ourselves here to a few very 
general remarks. The most interesting cases are then discussed in detail in 
Sect. |. 

The existence of l/Q power corrections in observables related to on-shell 
Green functions can already be seen from the simplest example, the two-point 
function of a heavy quark field, E(p, to). The 1-loop contribution is, schemat- 
ically, 

f d i k N 
E(p, to) cx J (27r)4 k2{2p , k + k2 + [p2 _ m2]) • (32) 

Only the denominator is important for the following discussion. Off-shell (p 2 ^ 
to 2 ) the contribution to the integral from k < A is of order CA 2 /to 2 relative 
to the dominant contribution from k ~ to. As p 2 — > m 2 the coefficient C 
diverges. Indeed, at p 2 = to 2 the integrand of Eq. (|32] ) vanishes only linearly 
with k for small k; hence we expect a non-perturbative contribution to E(m, to) 
of relative order A/m. 
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Similar power counting arguments apply to hadronic event shape observ- 
ables. The basic parton emission processes in QCD are infrared divergent for 
soft and collinear emissions. Event shape observables are constructed to sup- 
press these soft and collinear contributions, so that they are infrared finite 
in perturbation theory. But the residual contribution from small momentum 
regions is usually suppressed only by a single power of the hard scale. In gen- 
eral high-energy processes involving massless quarks the infrared contributions 
can be classified as soft or hard-collincar. It appears, however, that power 
corrections from hard-collinear regions (energy lu much larger than transverse 
momentum k±) are always suppressed by powers of Q 2 rather than Q. We do 
not know of a proof of this statement, but the following heuristic argument 
may illustrate the point: let p be the momentum of a fast on-shell particle, 
p ~ Q, after emission of a hard-collinear on-shell particle with w <~ Q, k± <C Q 1 
where k± is the transverse momentum relative to p. Then the propagator 



is expanded in fc^/w 2 ~ kj_/Q 2 and Q enters only quadratically. Since the same 
is true of the hard-collinear phase space, it may be argued that the transverse 
momentum, and hence Q, always enters quadratically as long as energies are 
large. As a consequence, if 1/Q power corrections exist and if one is interested 
only in those, the analysis simplifies, because only soft contributions need to 
be considered. 

A systematic analysis of IR renormalons for general short-distance pro- 
cesses has therefore much in common with the analysis of IR finiteness with 
the methods of perturbative factorization. It extends the notion of IR safety 
(absence of logarithmic divergences) to that of IR sensitivity (power-suppressed 
IR contributions). From a very general point of view, the most important les- 
son drawn from IR renormalons is the existence of a correlation between the 
size of non-perturbative corrections and the size of perturbative coefficients in 
large orders, often already at 2 loops. 

2.5 The Landau Pole 

In closing this overview section we address a common misunderstanding that 
the existence of renormalons and summation ambiguities is related to an in- 
frared Landau pole in the running coupling. A consequence of this misunder- 
standing is that it is often thought that the power corrections identified via IR 
renormalons have something to do with the definition of the QCD coupling. 
It is true that the factorial divergence follows from the fact that the coupling 
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evolves (/3o 7^ 0). The power correction indicated by this divergence is, how- 
ever, a property of the particular observable and it would exist even if the 
coupling did not evolve. 

To see how the association with the Landau pole might arise, we inter- 
change the summation and integration in Eq. (^5|). The interchange is math- 
ematically not legitimate, but if we nonetheless proceed, we obtain 
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where k 2 = —IQ 2 and 
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F{l)a s (ke' b ' 6 ), (34) 



as{ 1 l + a s ( M )/3oln(fc 2 /M 2 ) " falntf/A 2 { ' 

is the 1-loop running coupling. The Landau pole of the coupling lies on the 
^-integration contour and the ambiguity in defining the integral ([54]) due to 
this Landau pole is exactly identical to the ambiguity in defining the Borel 
integral of the divergent series expansion. 

This apparent connection does not persist beyond the calculation of bub- 
ble diagrams. The reason is that in the same approximation in which we 
keep only the bubble diagrams, the /3-function can only have a single term: 
/3(a s ) = —/3oa 2 . The existence of a Landau pole is then an automatic conse- 
quence. The general theory of renormalons shows that the leading asymptotic 
behavior depends only on (3q and fa, see Eqs. ( p^|Jl7| ). On the other hand 
whether a Landau pole exists or not is a strong-coupling problem and it de- 
pends on all coefficients of the /3-function, and on power corrections to the 
running of the coupling. This pari hp studied explicitly for /3-functions with 
only two independent terms.E3&L-3 For example, the integral ( pi) ) might be 
well-defined, but its series expansion remains still divergent. It would then be 
incorrect to conclude that the power corrections indicated by IR renormalons 
are properly taken into account by summing the series to the numerical value 
given by the integral. Such a summation prescription would be related to 
a certain definition of non-perturbative parameters relevant to the particular 
process, but it would not render these parameters zero and superfluous. In 
essence, IR renormalons reflect perturbative aspects of non-perturbative cor- 
rections (namely their power divergences) for a specific observable, whereas the 
existence of a Landau pole is a wholly non-perturbative issue, but not related 
to any particular observable. This can again be non-perturbatively verified in 
the non-linear O(N) <7-model,t£l where the propagator of the a-field defines an 
effective coupling without Landau pole, but all correlation functions contain 
the full series of IR renormalon poles. 
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3 Applications 

3.1 Deep Inelastic Scattering 

We begin with a short exposition of the applications of renormalons to the 
structure functions of deep inelastic scattering. Since the operator product 
expansion is available, this example serves to illustrate the operator interpre- 
tation of renormalons and formulate phenomenological procedures that can be 
generalized to other, more complicated situations. 



The GLS Sum Rule 

The GLS sum rule, defined in Eq. ([j]), has already been discussed in the intro- 
duction; here we complete this discussion. 

The first IR_renormalon singularity at t — 1//3q corresponds to a single 
twist-4 operators and can be obtained with the methods described in Sect. |[ 
Combining this result with the calculation of the leading ultraviolet renor- 
malonE£l (see Sect. ||), we find the large-order behavior 



ifUV (_ 1)nn l-fSi/ft+\, + R m n /3 1 //3 2 -32 b /9l i (36) 



where r n is the coefficient of /?o,i are the first two coefficients of the j3- 

function, b — 47r/?o = 11 — 2iV//3 and Ai is related to the anomalous dimension 
matrix of four-fermion operators, see Sect. 4. For Nf > 2, the UV renormalon 
behavior dominates the asymptotic behavior at very large n because of its 
larger power of n. However, the overall normalization constants K uv and K m 
are not known. Since the MS scheme favors large residues of IR renormalons, 
at least in the large-/?o approximation, we expect fixed-sign IR renormalon 
behavior in intermediate orders. The first three terms in the series known 
exactly are indeed of the same sign in the MS scheme as can be seen from 
Eq. #• 

Is the asymptotic behavior in Eq. ( |36| ) r eleva nt to phenomenology? Since 
the constants K are not known, we consideiOtZl the large-/?o approximation. 
The Borel transform (defined by Eq. (^)) of the perturbative expansion in this 
approximation is given byEa 

where u = fat. This is much simpler than what we would have expected on 
general grounds. In particular, there are only four renormalon poles, all others 
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being suppressed in the large- (3$ limit. The structure of the leading singular- 
ities is also simpler than the exact result in Eq. (|3^), because the anomalous 
dimensions should be set to zero in this limit. Nonetheless, the exact n, ri are 
reproduced reasonably well by the large-/?o limit. The large-/?o approximation 
taken at face value implies that the minimal term of the series is reached at 
order or af at Q 2 = 3GeV 2 , a momentum transfer relevant to the CCFR 
experiment. Hence it is not clear whether at Q 2 — 3 GeV 2 the perturba- 
tive prediction could be improved further by exact calculations of higher-order 
corrections. Further improvement would then require the inclusion of twist- 
4 contributions, and in particular a practically realizable procedure to com- 
bine them consistently with the perturbative series. This hypothesis is further 
supported by noting that the integral over loop momentum is dominated by 
k ~ 450 MeV at order and k ~ 330 MeV at order af. The ambiguity in 
summing the perturbative expansion is of order (assuming Aj^t =215 MeV) 

1 „ 1 8e 5 / 3 A§^ 0.1 GeV 2 

g SiGLS - J^~Q^~ ~Q2 Q2 ■ ( 38 ) 

This should be compared to the twist-4 contribution to the same quantity 
estimated in quark model and by QCD sum rulesLj 

8 ((Q 4 )) _ 0.1 GeV 2 

27 Q 2 ~ Q 2 ' [ ' 

where ((O4)) is the reduced nucleon matrix element of the twist-4 operator. 
The two are comparable, which suggests that the treatment of perturbative 
corrections beyond those known exactly is as important for a determination of 
a s from the GLS sum rule as the twist-4 correction. 



The 2>Dependence of Power Corrections to Structure Functions 

The operator product expansion allows us to express 1 /Q 2 corrections to the 
structure functions in terms of contributions of several towers of twist-4- oper- 
atorsH or, equivalently, several multiparton correlation functions. circd'O Un- 
fortunately, the structure of the corrections is complicated and they involve 
many non-perturbative parameters which cannot all be extracted from inclu- 
sive measurements. Because of this, it has never been possible to use this 
sophisticated machinery in the analysis of real data. In practice, higher-twist 
corrections are being extracted from data from a combined fit to perturbative 
and l/Q 2 -suppressed contributions in a large Q 2 range. 

Renormalons provide a simple ansatz for the cc-dependence of power cor- 
rections involving fewer (if any) non-perturbative parameters. As in most 
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other applications divergences of the perturbative series are relevant inasmuch 
as they originate from small momentum regions in Feynman diagrams that are 
also responsible for non-perturbative effects. To explain the idea, recall that 
the structure functions (take ^2(2:) as an example) are related to the parton 
distribution functions P(x), P = q,q,g through a perturbative expansion and 
convolution of the form 



F 2 (x 7 Q 2 ) 



2x 



+ 



p 



el 



n=0 



r PiV ,(x)a 



n+l 



cp{x)- 2 



*P(x,Q 2 ) 



C 2 (x) 

Q 2 



(40) 



where ep are the electromagnetic charges of partons and * stands for the 
convolution. In this expression we have subtracted contributions of small mo- 
menta k < [i <C Q from the coefficient functions and defined the higher-twist 
contribution which can have a complicated operator content as the full con- 
tribution (perturbative and non-perturbative) coming from small momenta 
k < ji. Since twist-4 operators are quadratically divergent, we know that 
C 2 (x) = c'(x)^ 2 + 0(A 2 ) if > A and c'(x) = J2p cp(x) * P(x) so that the 
dependence on the cutoff /i cancels. 

As c(x) is defined as the small-momentum contribution to Feynman dia- 
grams, it has itself a (divergent) perturbative expansion in a s . We then extract 
the c(x) from the x-dependence of the IR renormalon pole computed in the 
large-/?o approximation. In this .appr.oximation there is only a (anti-)quark 
contribution and the result isEScaelJ'C^I 



cf\x) 
cf\x) 



S\x) 



8x 2 - 45(1 -x), 



(41) 



[1-x], 
4 

[1 - A- 



+ 4 + 2x+ 12x 2 -9<5(l-x) -S'(l-x), (42) 
+ A + 2x + Ax 2 -56(l-x)-6'(l-x) (43) 



for Fl, F 2 and F3, respectively. A common overall normalization is omitted 
here, because it plays no role in the following. The '+' prescription is defined 
as usual by J„ dx [f(x)]+t(x) = J Q dx f(x) (t(x) — t(l)) for test functions t(x). 

If we assume that the subtraction term in the square brackets of Eq. (40) 
cancels approximately the higher-order perturbative terms, and if we assume 
that C 2 (x) is approximated by c'(x)[i 2 , we obtain an improved prediction for 
the structure functions compared to the fixed-order perturbative approxima- 
tion. This suggestion, though sometimes motivated by different arguments, has 
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become known as "renormalon model" of power corrections .c3Ej The structure 
functions are then written as 

F(x,Q) = F*"- 2 {x,Q) (l + + 0(1/Q 4 )) , (44) 

where F tw ~ 2 (x, Q) is the leading-twist result for the structure function F — 
F L ,F 2 ,F 3 ,... and 

A 2 r 1 d£ 

D(X > Q) = F*>-*{x,Q) I J C(0 q{x/ ^ M) (45) 

is the model parametrization of the (relative) twist-4 correction. Here q(x, //) 
is the standard (leading-twist) quark density, and A is a certain scale of order 
Aqcd which provides the overall normalization. The expression can be ex- 
tended to include gluon contributions and/or higher-order corrections should 
they become available. 

The overall normalization has been treated differently in the literature. 
One suggestion has beenc3 to parametrize the normalization of all 1/Q 2 power 
corrections by a single process-independent number, to be extracted from the 
data once and related to a a certain effective QCD coupling. Other au- 
thorsOE-jO prefer to adjust the normalization in a process-dependent way 
and to take only the shape of the x-distribution as a prediction of the model. 
Because of difficulties in constructing the gluon contribution in the model, one 
may also think of adjusting the normalization of quark and gluon contributions 
separately. 

The "renormalon modal!! .of. twist-4 corrections has first been applied to 
the structure function F2.CJO'c£l As shown in Fig. ||, the shape of the twist-4 
correction calculated from the model indeed reproduces the experimental data 
very well. These results refer to the the non-singlet contribution to F2, which 
is expected to dominate except for small values of x. Similar p redictio n s have 
been worked out for the longitudinal structure function i^,E2la i^3,e3o and 
the polarized structure function g\ etc.r 9 ^ 1 ! More recently a renornialoa model 
prediction has also been constructed for the singlet contributiorBS to F2, 
which modifies the analysis at small x, below those x for which comparison 
with present data is possible. (The treatment of singlet contributions is more 
difficult and ambiguous in the renormalon model than non-singlet contribu- 
tions. The calculation relies on singlet quark contributions, which are then 
reinterpreted as gluon contributions.^!!) 

A striking property of the renormalon model for twist-4 corrections is 
that all target-dependence enters trivially through the target dependence of 
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Figure 3: Relative twist-4 contribution D(x) (called C Pt d{x) here) denned WiEq. ( fli) ) 
to the proton (deuteron) structure function Fi in the "renormalon model"c3 (dashed 
line) compared with proton (filled circles) and deuteron (empty circles) data.E-3 The 
solid curve shows the unrescaled estimate of the renormalon ambiguity. 



the twist-2 distribution functions. Therefore the renormalon model can be 
useful only if the genuine twist-4 target dependence is small compared to the 
magnitude of the twist-4 correction itself. In terms of moments M n , Eq. (45) 
implies 

/M-tw-4 M tw ~ 4 



Ihadronl 



Mn 



|hadron2 



(46) 



,ust deuterons, 
3 also confirm 



Figure || shows that this is indeed the case for F2 of protons ag; 
in particular in the region of large x. More recent analysescS 
this fact. 

It is known that higher-twist corrections (as well as higher-order perturba- 
tive corrections) are enhanced as x — > 1£3 This is in part an effect of kinematic 
restrictions near the exclusive region and the renormalon model reproduces 
such enhancements. For the structure functions it is found that power correc- 
tions related to renormalons are of order 



A 2 



(47) 



f(l-x)_ 

at least those related to diagrams with a single gluon linc@ This tells us that 
the increase of the twist-4 correction towards larger x seen in the model and 
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-1 '—J ' — ' — 1 

-2 -1 
10 10 X 1 

Figure 4: Twist-4 correction to xF$ as extracted from the (revised) CCFR data. The 
three plots show the effect of including leading order (LO), next-to-leading order 
(NLO) and next-to-na*jt-to-leading order (NNLO) QCD corrections in the twist-2 
term. The data pointsEZI are overlaid with the shape obtained from the "renormalon 
model" for the 1/Q 2 power correction. 
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the data in Fig. |^ may to a large extent be the correct parametrization of such 
a kinematic effect. Note that Eq. (|47]) can be understood from the fact that 
the hard scale in DIS is — x at large (but not too large) x. 

It is also possible that both the experimental parametrization of higher- 
twist corrections and the model provide effectively a parametrization of higher- 
order perturbative corrections to twist- 2 coefficient functions. As far as data 
are concerned, it should be kept in mind that it is obtained from subtracting 
from the measurement a twist-2 contribution obtained from a truncated per- 
turbative expansion. As far as the renormalon modal is concerned, it is best 
justified by the "ultraviolet dominance hypothesis" O Since UV contributions 
to twist-4 contributions can also be interpreted as contributions to twist-2 
coefficient functions, a "perturbative" interpretation of the model prediction 
suggests itself as we have already indicated in the discussion of Eq. (p0|) . Note 
that higher-order corrections in a s (Q) vary more rapidly with Q than lower 
order ones, and may not be easily distinguished from a 1/Q 2 behavior, if the 
Q 2 -coverage of the data is pot-spry large. An interesting hint in this direction 
is provided by the analysisE3EK*l of CCFR data on F3, reproduced in Fig. |[ 
The figure shows how the experimentally fitted twist-4 correction gradually 
disappears as NLO and NNLO perturbative corrections to the twist-2 coef- 
ficient functions are included. At the same time, the renormalon model for 
the twist-4 corrections reproduces well the shape of data at leading order, and 
hence parametrizes successfully the effect of NLO and (approximate) NNLO 
corrections. This is an important piece of information, relevant to quantities 
for which an NNLO or even NLO analysis is not yet available. 

Note that whether the model is interpreted as a model for twist-4 cor- 
rections or higher-order perturbative corrections is insignificant inasmuch as 
renormalons are precisely related to the fact that the two cannot be separated 
unambiguously. The model clearly cannot be expected to reproduce fine struc- 
tures of twist-4 corrections. Its appeal draws from the fact that it provides a 
simple way to incorporate some contributions beyond LO or NLO in perturba- 
tion theory, which may be the dominant source of discrepancy with data with 
the presently achievable accuracy. 

3.2 Hadronic Event Shape Variables 

The structure of hadronic final states in e + e~ annihilation and deep inelas- 
tic scattering is the subject of intensive ongoing studies. E3 This structure is 
characterized by a set of infrared and collinear safe event shape variables that 
are calculated in perturbative QCD in terms of quark and gluon momenta and 
compared to the measured hadron distributions. Apart from a correction for 
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detector effects, the comparison of theory and data therefore requires a correc- 
tion for hadronization effects that are most commonly modeled using Monte 
Carlo event generators. It has been known for quite some timeEj that the 
hadronization corrections are substantial. In this section we review recent de- 
velopments that relate hadronization corrections to power corrections of order 
A/Q (where Q is the center-of-mass energy in e + e _ annihilation) indicated by 
renormalons in the perturhative prediction for the event shape variables. This 
connection was suggestedEE3 and worked out for a few cases of practical inter- 
est!^ several years ago. These studies provided the first theoretical indications 
that hadronization corrections to most of the observables should scale as A/Q, 
i.e. are surrojjessed by only a single power of the large momentum. Subsequent 
analyses! 6 •ffi 3 rEjE3 confirmed this conclusion. 

Most of the discussion below assumes a generic event shape observable S 
that is of order a s (Q) at leading order. One example is S = 1 — T where 
"thrust" T is defined as 

r = maxE^, (48) 

Li \Pi\ 

where the sum is over all hadrons (partons) in the event. The thrust axis Ht is 
the direction at which the maximum is attained. Other event shapes considered 
in connection with power corrections are the heavy jet mass, jet broadening, 
C-paramctcr and the longitudinal cross section in e + e~ annihilation, to give 
only a few examples. 

Mean Values of Event Shape Variables 

It is relatively easy to understand that event shape observables in e + e~ annihi- 
lation arc linearly sensitive to small parton momenta and hence are expected to 
receive non-perturbative contributions of order A/Q. Consider an event shape 
variable that is zero at tree level and therefore related to the matrix element 
for gluon emission 7* — ► qqg to leading order 

(S) = JdPS[ Pl }\M q g g \ 2 S( Pt ). (49) 

It can be argued that the A/Q sensitivity arises neither from emission of 
collinear and energetic partons, nor from soft quarks, but only from soft glu- 
ons. Introducing the energy fractions Xi — 2pi ■ q/q 2 , and reserving 23 for the 
gluon energy fraction, this implies that the only relevant integration region 
is 23 — > and, therefore, the gluon emission can be calculated in the "soft" 
approximation: 

^ 2 = 329 ° (i- Xl Ui- X2 y (50) 
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The phase space is, on the other hand 



dPS[pi] = / dx\dx2 0(xi + x% — 1) 



(51) 



and since the matrix element in Eq. fl50| ) is singular as — > 1, there is a 
potential logarithmic singularity. To obtain an IR finite result, the event shape 
variable has to be constructed so as to eliminate this divergence. The generic 
situation with event shapes is a linear suppression of soft gluons, S{xi) oc X3 
as X3 — > 0, e.g. for the thrust 1 — T = 1 — max(xi, X2) with x\ + X2 + £3 = 2. 
It is easy to see that this property implies a contribution of order fi/Q to the 
integral in Eq. (^) from gluons with energy less than fi, unless there is some 

cancelation. |—||—| r~ 1 

It has been suggestedE3S'E3 that the leading power correction to average 
event shape observables may be described by a single ( "universal" ) parameter 
multiplied by an observable-dependent, but calculable, coefficient. Write 



(S) = A s a s (v) + 



2 1 



B s + A s Po In ^ 



Q 



0(1/Q 2 ). (52) 



Dokshitzer and Webber 
in the form 



parametrize the coefficient of the power correction 
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where fij is an IR subtraction scale (typically chosen to be 2GeV), ao (/•*/) is 
the non-perturbative parameter to be fitted and K = 67/6 — 7r 2 /2 — 5JV//9. 
The remaining terms approximately subtract the IR contributions contained 
in the perturbativc coefficients A and B up to second order. The universal- 
ity assumption can be tested by fitting the value of ao(/Ur) or i equivalcntly, 
Ks{[i)/cs to different event shape variables. 

In Fig. ^| we compare the energy dependence of (1 — T) and the heavy jet 
mass (Mfi/Q 2 ) with the prediction with and without a 1/Q power correction. 
It is clearly seen that (a) the second-order perturbative result with scale fi = Q 
is far too small and (b) the difference with the data points is fitted well by a 
l/Q power correction. 

In absolute terms the power correction added to thrust and the heavy jet 
mass is about 1 GeV/Q. This is a sizeable correction of order 20% even at the 
scale Mz, because the perturbative contribution is of order a s {Mz)/'K. The 
fit for cio is sensitive to the choice of renormalization scale fi and in general 
to the treatment of higher-order perturbative corrections. There is nothing 
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Figure 5: Energy dependence of (1— T) (upper panel) and the heavy jet mass (Mjj /Q 2 ) 
(lower panel) plotted as function of l/Qx-\ Dotted line: second-order perturbation 
theory with scale \i = Q. Solid line: second-order perturbation theory with power 
correction added according to Eq. j53| ) and with = Q, fii = 2 GrpV. For d?o(2GeV) 
the fit values 0.543 for thrust and 0.457 for the heavy jet masa-J are taken. (Note 
that this reference uses Ci_t = c M 2 i s = 1.) The dashed line shows second order 
perturbation theory at the very low scale 0.07Q with no power correction added. For 
both observables a a (Mz) has been fixed to 0.12. 
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wrong with this, because the very spirit of the renormalon approach is that 
perturbative corrections and non-perturbative hadronization corrections are to 
some extent inseparable. Hence we find it plausible that the 1/Q power cor- 
rection accounts in part for large higher-order perturbative corrections, which 
are large precisely because they receive large contributions from IR regions of 
parton momenta. It has been notedEa that choosing a small scale, [i = 0.13Q, 
reduces the second order perturbative contribution and power correction sig- 
nificantly for (1 — T). In Fig. || (dashed curve) we have taken a very low scale, 
fi = 0.07Q, to illustrate the fact that the running of the coupling at this low 
scale can fake a 1/Q correction rather precisely (a straight line. in the figure). 
Note that an analysis of (1 — T) in the effective-jcharge schemeL^I selects almost 
the same scale fj, — 0.08Q. A simultaneous fit£3 of a s , a third-order pertur- 
bative coefficient and a 1/Q power correction then leads to reduced power 
correction of order (0.3 ± 0.1) GeV/Q consistent with the above argument. 

A comprehensive compilation^] of the combined fits of the the parameter 
ao(fj,i) and the strong coupling as (Afz) to the data on several event shapes 
in e + e~-annihilation and deep inelastic scattering is shown in Fig. |^. The 
extracted values of ao center around <5o(2 GeV) ~ 0.5 that gives some sup- 
port to the universality hypothesis. Its theoretical status has not been com- 
pletely elucidated so far, different and somewhat conflicting arguments have 



30 



been givenE3SE3 In order to get further insight, a detailed analysis of IR 
sensitive, contributions to the matrix elements for the emission of two par- 
tonsO'Ol has been undertaken. It was found that for (1 — T), the jet masses, 
and the C-parameter the coefficient of the 1/Q power correction that is ob- 
tained for one gluon emission is rescaled by the same factor, often called the 
"Milan factor", whose value was later revised to 1.4900 The universality 
of the Milan factor for a class of event shapes is a direct consequence of the 
assumed dominance of soft-gluon radiation coupled with an underlying geo- 
metrical universality (linearity in transverse momenta of emitted gluons) in 
the shape variables themselves.L3 Note that the numerical value of the "Mi- 
lan factor" is actually not important, since it is likely to be modified by yet 
higher-order corrections. The universality of the "Milan factor" is hence the 
more important observation, since it allows us to relate different observables. 
Nevertheless, the analysis of soft-gluon effects at the 2-loop order is very in- 
teresting, since it is represents one of the few cases where power corrections 
have been explicitly investigated beyond one-gluon exchange for observables 
that do not have an operator product expansion. 



Event Shape Distributions 

The structure of power corrections to event shape distributions (l/o~tot)do-/dS 
is more complex. (The distributions are defined so that the mean value is 

given by (S) = (1/crtot) J man: dS S da / dS .) In the following we consider event 
shapes defined so that S ^ corresponds to the two-jet limit. For small values 
of the event shape variable S the dynamics of soft-gluon emission depends on 
two different IR scales QS and y/Q^S, and 1/Q < l/Qy/S < 1/(QS). The 
smallest scale of order QS is related to the typical energy carried by soft gluons, 
while the scale Qy/S defines the transverse momenta of the jets, k\ ~ Q 2 S. By 
examining the sensitivity of perturbative emission of soft gluons with energy 
of order QS and collinear particles with transverse momentum of order Q 2 S, 
we are lead to suspect non-perturbative corrections suppressed by powers of 
both scales. Then, since in the end-point region S ~ A/Q, we can expand 
the distributions in powers of the larger scale Q 2 S and keep the leading term 
only that corresponds to the resummation of all corrections of order (l/QS) k , 
neglecting corrections of order (l/Q\fS) k . The resummation introduces the 
important concept of a shape function for soft-gluon emission,L3'LD which we 
shall review briefly on the particular example of the thrust distribution t = 

i-tE 

The starting point is the observation that the differential thrust distribu- 
tion da/dt in the small t region computed by resumming an infinite number of 
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soft-gluon emissions is expected to exponentiate under a Laplace transforma- 
tion 

— f^dt e -» t ^ = e- s ^ (54) 

&tot Jo at 

with the exponent S(v, Q) of the general form 

S(y,Q) = /^(l-^) 
Jo x 

x / -^T[a s (kl)} + B[a s (xQ 2 )} + C[a s (x 2 Q 2 )} . (55) 

\Jx 2 Q 2 k ± J 

Here T[a s ] is the universal cusp anomalous dimension that controls soft and 
collinear gluon emission. (Compared to the original discussion^ we have 
added the additional function C[a s (x 2 Q 2 )] to the exponent. This function 
does not appear to the next-to-leading logarithmic accuracy, but we are not 
aware of a theoretical argument that would not allow terms of this structure 
in higher-orders. The following discussion proceeds under the assumption that 
C[a 8 (x 2 Q 2 )] does not have a ^divergent perturbative expansion, as it can hap- 
pen in Drell-Yan productiono) The next step separates the contribution of 
soft gluons introducing a cutoff in transverse momentum k± < fi. Heuris- 
tically, the contribution of gluons with k± > fi has to be defined as a per- 
turbative contribution Spt(v,Q, fj) to S(v,Q), and contributions of k± < fi 
have to be promoted to the non-perturbative correction, Snp(v, Q, /•*), so that 
S{v,Q) = Spt(v,Q,ij)+Snp(v,Q,h). From the T-term in Eq. (p5|), expanding 
formally in powers of v /Q and neglecting all powers of v/Q 2 , we obtainLj 

n ^ 

SNpfaQ,?)"^—^^) I dklk n - 2 Y[a s (kl)]. (56) 
^ rml V Q J J 



(i 



The integrals over the cusp anomalous dimension should eventually be substi- 
tuted by (dimcnsionful) non-perturbative parameters. Since the variable v is 
conjugate to t, Eq. ( |56|) effectively organizes all power corrections in l/(tQ). 

If t 3> A/Q, then keeping the first term only in the sum in Eq . (p6[) is 
sufficient, S^p^, Q, fi) ~ \\v /Q. Assuming exponentiation as in Eq. (|54j), the 
non-perturbative correction ampunts in this case to a shift in the resummed 
perturbative thrust distributior£§Ej 

% (t) — > ^(t - Xi/Q) + 0(l/(tQf) , (57) 
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with Ai the same non-perturbative parameter that parametrizes the 1 /Q cor- 
rection to the mean thrust, i.e. (i)i/<g ~ Ai/Q. 

If, on the other hand, t ~ A/Q, then all terms in the sum in Eq. (|56| ) have 
to be kept. The infinite set of non-perturbative parameters corresponding to 
the increasing powers of v /Q defines a Q-independent function /t(e; fi) through 

/■oo 

/ de exp(-z/e/Q)/ t (e;/i) = eyvp(-S N p{Q/v, //)). (58) 
Jo 

The function / t (e;/x) parametrizes the energy spectrum of non-perturbative 
soft-gluon emission for the thrust, observable, similar to the shape function 
introduced for the descriptionP°tP 9 rr°l of the endpoint region in inclusive heavy 
meson decays such as B — > jX s . With this definition, the differential thrust 
distribution corrected for the non-perturbative effects takes the forrr£3 

— ^r = QMtQ)Rp T (0,fi)+ [ tQ def t (e^) — ^(t-e/Q), (59) 
otot at Jo otot at 

where R(0, /i) is the Sudakov factor taking into account the contribution of 
virtual soft gluons with the energy above the cutoff /i and the subscripts "PT" 
indicate that the corresponding quantity is calculated in perturbation theory. 
The shape function induces a smearing of perturbative gluon radiation by non- 
perturbative corrections. Note that shape functions depend, in general, on the 
particular event shape observable. In addition to thrust, &\so the heavy jet 
mass and the C-parameter distributions have been studiedO 

In contrast to heavy quark decay where non-perturbative corrections ex- 
tend the photon spectrum beyond the perturbative boundary of phase space 
due to the energy distribution of the heavy quark in the meson, the non- 
perturbative corrections to the thrust distribution have an opposite effect. 
They shift the distribution inside the perturbative window < t < t max and 
describe the "evaporation" of the energetic jets in the final state, loosing en- 
ergy for the emission of soft gluons. In practice, one has to model the shape 
function using a certain ansatz and fit the parameters to describe the thrust 
distribution at a certain value of the cm. energy s = Q 2 . The result of such a 
fit is shown in Fig. [?]. The energy dependence of the differential thrust distri- 
bution is then predicted without free parameters and appears to be in a good 
agreement with all available data, see Fig. || 

The Energy Flow Correlation Function 

The most important statement that has emerged from the application of renor- 
malons to event shapes and is supported by all existing evidence is that soft- 
gluon emission presents the only source of linear A/Q power corrections. This 
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Figure 7: Fitting the shape functionB from the comparison with the data on the 
differential thrust distribution at Q = 91.2 GeV. 
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Figure 8: Using the shape functiorH to predict the differential thrust distribution at 
Q = 14, 22, 35, 44, 55, 91, 133, 161 GeV (from bottom to top). 
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observation has profound consequences since emission of soft gluons occurs at 
time scales that are much larger than those involved in the formation of per- 
turbative narrow jets. As a consequence, these two subprocesses are quantum- 
mechanically incoherent and the soft-gluon distribution emitted by a pair of 
quark jets at wide angles depends only on the direction and total color charge 
of the jets. This heuristic reasoning suggests factorization: to all orders in 
perturbation theory, inclusive cross sections for two narrow jets in e + e~ anni- 
hilation can be written as products of separate functions for the jets and the 
soft-gluon radiation, up to corrections suppressed by powers of (A/Q) 2 . As 
a consequence, we can use the eikonal approximation and replace the jets by 
eikonal, or Wilson lines W± — Pexp[ig s L dup±A tl {up±)) in which soft-gluon 
field is integrated along the light-like directions p± defined by the momenta of 
the outgoing jets. In this approximation the shape function for a given event 
shape variable S becomes 



where the sum goes over all final states and an UV cutoff \x is assumed for the 
gluon momenta. 

The sum over the final states can be performed by introducing an operator 
£{ft) that measures the_density of energy flow in the direction of the unit vector 
ft at spacial infinityEpp Then, in particular, the shape function for the thrust 
distribution becomesEil 



f t {e) = {0\{W + Wl){0)5[e- dnw t {ft) 8 {ft, p) (W-Wj) (0)|0), (61) 



where uit {ft) contains the information on the particular shape variable. In gen- 
eral, the complete information about soft-gluon emission is encoded in "multi- 
energy flow" correlation functions 



g{ft 1) ... ) ft N ; f i) = {0\{W + wl){Q)£{n 1 ),...,S{fi N ){W-Wl){0)\Q) (62) 



that describe the energy flow at spatial infinity. Power corrections to different 
event shape averages can be calculated in terms of the energy flow as 



The concept of energy flow is very useful, since it allows us to define the quan- 
tities of interest on an operator level, which makes them more amenable to a 




(60) 





(63) 
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rigorous analysis of power corrections. On the other hand, many important 
issues remain to be resolved such as the dependence of the energy flow correla- 
tion functions on the factorization scale fj,. But already at the present stage the 
concept of energy flow has provided insights into the extent to which universal- 
ity can be expected for power corrections to various moments or distributions 
of event shape observables. 

3.3 Heavy Quarks 

In this section we consider hard processes for which the large scale is given 
by the mass of a heavy quark. We discuss the notion of the (pole) mass of a 
heavy quark itself and its relation to the heavy quark potential, and consider 
applications of these results to tt pair production near the production threshold 
in e + e~ annihilation. 

The Pole Mass 

Since quarks are not observed as asymptotic states, their masses generally have 
to be considered as parameters in the QCD Lagrangian, on par with the strong 
coupling. Because the running of the coupling is conventionally considered 
using dimensional regularization, it is most natural to also employ the MS 
scheme for the mass definition, introducing running quark masses fh q (fi-^) 
and fixing their values at a certain reference scale. This is indeed the procedure 
used to deal with light quarks q = u,d,s and also heavy quarks provided the 
hard scale in the process is larger than or of order of the quark mass. On the 
other hand, using MS heavy quark masses is not convenient in processes where 
the hard scale is significantly smaller than the mass of the quark itself. The 
reason for this is that the usual renormalization group expression 

\a s (m b )J 

is physically irrelevant at fi <SC nib because it is derived by assuming that /j. is 
the UV cutoff and thus the largest scale. (There is, formally, nothing wrong 
with taking /i <C nit, in Eq. (|64|). However, in calculations of physical observ- 
ables, e.g. heavy quark decay rates, we expect large perturbative corrections in 
higher orders in this case, because unphysical, large logarithms are generated.) 
On the other hand, in the limit fi <C mj, the heavy quark interacts with gluons 
through the color Coulomb potential V(r) = — Cfa s (l/r)/\f\. A restriction 
on the gluon momentum \k\ < (i corresponds to the cutoff at large distances 
r > 1/fJ, so that the, dependence of the mass parameter on /i at A <C [i <C m\, 
is in fact linearE3E3H 
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Figure 9: A schematic representation of the scale dependence of a heavy quark mass 
defined with a "physical" IR cutoff (solid line) and in the dimensional regularization 
(dashed line). The shaded area illustrates uncertainties of the mass definition when 
the scale is driven to values of order A 



m b (0) - m b (fi) — Cfa s n , 



(65) 



see Fig. ^| for an illustration. The quantity m(0) would correspond to a "phys- 
ical" quark mass if it existed and in perturbation theory can be defined as the 
location of the pole in the perturbative quark propagator, i.e. the pole mass: 



"Ipole — TO Ms( m Ms) 



1 + H r " a "( TO Ms) 



(66) 



The pole mass is I 
mali jajjkm^cheme £ 



iriite, gauge independent and independent on the renor- 
3 The perturbative series in Eq. (pq) is, however, diver- 



r n ~ const • (2#,) n r[n + 1 + A/(2/3 ) 5 



(67) 



The sum of the series is ambiguous by an amount of order A, and, therefore, 
the quark pole mass is perturbatively defined only to an accuracy 



#TOpole ~ A. 



(68) 



This uncertainty has important practical implications as it means that the 
pole mass has to be eliminated as a parameter in calculations of physical ob- 
servables, if these observables are less sensitive to the IR region than the pole 
mass itself. Important examples of such observables are inclusive heavy quark 
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decaysBBElIliJEil and top quark production near threshold in e + e annihi- 
lation. The second example will be discussed in some detail below. IR_ceiior- 
malons have also been investigated for exclusive heavy quark decaysJSEj in 
which a relation with the heavy quark pole mass emerges through the binding 
energy A, defined in heavy quark effective theory. 



The Heavy Quark Mass with IR Subtractions. 

In order to deal with a well-defined non-perturbative parameter and at the 
same time to retain the attractive features of the pole mass we have to per- 
form an explicit scale separation. Hcuristically, this would mean that the 



perturbative contributions in Eq. d66|) have to be calculated using a cutoff /i 
at small momenta. The dependence of such subtracted mass on fi is linear and 
the coefficient depends on the particular procedure to implement the IR-cutoff, 
but it is difficult in practice to implement such a procedure, and to maintain 
gauge invariance in_particular. 



One suggestion^ has been to utilize a close connection between ambiguities 



in t 
tior 



S3 



pole mass and the static heavy quark potential. The starting observa- 
is that the leading IR power correction to the potential in momentum 
space cannot be A/|<f|, but has to be quadratic: 

nt) = - 4?rC y (<f) (l + • ■ . + const •!? + •••)■ (69) 

To avoid misunderstanding, note that we are not concerned with the long- 
distance behaviour of the potential at q ~ A, but with the leading power cor- 
rections of the form (A/q) k , which correct the perturbative Coulomb potential 
when q is still large compared to A. 

When we consider the coordinate space potential, given by the Fourier 
transform of V(q), a new situation arises. It is easy to see by dimensional 
analysis that the contribution of small \q\ < [i in the Fourier integral is a 
r- independent constant: 

i9 ' r V{$)~ [ jLLytf) = const -fi. (70) 



(2tt) 3 J (2tt) 3 

l<fl<W l<?l<« 

This implies a long-distance and hence non-perturbative correction 

V(r) = - CFOls ( 1 / r ) (1 + , . . + const • Ar + ■ ■ .) (71) 
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for the coordinate space potential that can also be observed through the cal- 
culation of renormakms generated by one-gluon exchange with vacuum polar- 
ization insertions £3o 

The r-independent constant in the potential is closely related to the un- 
certainty in the pole mass. This is easy to understand since the total energy of 
a heavy quark-antiquark system E{r) = 2m po i c + T^Cd-is a physical observable 
and has to be well defined in the heavy quark limit Sa We can thexL. define a 
potential-subtracted (PS) quark mass and a subtracted potential asEJ 



where 



mpsifJ-f) = TOpole - Sm{p,f). 
V(f,n f ) = V(f) + 25m(n f ) 



\q \<Vf 



(72) 



(73) 



To 2-loop accuracy, the relation of the PS mass to m = Wj^g (m-^g) is given by 
_ J 4a s (m) 
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>, (74) 



where K\ = 13.44— 1.04A^ is the 2-loop coefficient in the relation of m po i to 
TTitB and K<2, — 10.33 — l.lXNf the 1-loop correction to the Coulomb potential 
in momentum space. (The 3-loop relation is also known.) The important point 
is that in the perturbative expansion that relates the two masses in Eq. (74) 
the leading IR renormalon divergence has been eliminated. In practice, this 
leads to smaller perturbative coefficients starting already at two loops. 

We can use the PS mass and subtracted potential instead of the pole mass 
and the Coulomb potential to perform Coulomb resummations for threshold 
problems. The benefit of using an unconventional mass definition is that large 
perturbative corrections related to strong renormalon divergence associated 
with the coordinate space potential are obviated. Physically, the crucial point 
is that, contrary to intuition, heavy quark cross sections near threshold are in 
fact less long-distance sensitive than the pole mass and the coordinate space 
potential. The cancelation is made explicit by using a less long-distance sensi- 
tive mass definition. 
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Top Quark Production in e + e Annihilation 

It is perhaps surprising that the renormalon divergence in the pole mass/MS 
mass relation and the mass subtractions discussed above have played a very 
important role for the top quark, which is so much heavier than the QCD 
scale. The reason for this is the extraordinary precision of about 100 MeV 
with which the top quark mass can in principle be determined by scanning the 
pair production threshold at an e + e~ collider. 

The improvement in convergence due to the subtraction term can be seen 
on the one hand in the perturbative conversion to the MS mass, and in the ti 
line shape on the other hand. Numerically, the series that relate the pole and 
PS mass, respectively, to the MS mass fh t = rh t (rh t ) are as follows: 

m t = [165.0 + 7.64 + 1.64 + 0.52 + 0.25 (est.)] GeV (75) 
m t , PS (20 GeV) = [165.0 + 6.72 + 1.21 + 0.29 + 0.08 (est.)] GeV. (76) 

for fht — 165 GeV and a s (rht) = 0.1091 (corresponding to a s ( m f)in 0.119). 
The 3-loop coefficients can also be comDuted using recent results. Ej-EJ The 4- 
loop estimate uses the "large- /V' limit Hj The improved convergence is evident 
and significant on the scale of 100 MeV set by the projected uncertainty on the 
mass measurement. 

The calculation of the top quark pair production cross section near thresh- 
old uses a non-relativistic effective theory to accomplish the all-order resum- 
mation of Coulomb "singularities" . We cannoL_go into details here except 
to mention that the leading order cross sectiorO can be obtained from the 
effective Lagrangian 

Cos = ^ ^d° + Sm t (M/) + J^+*^ ^ W 
+X t lid° + 5m t (fi f )-^—+iT t ) X 



rtyV](r)(-^) [X f x](0), (77) 



where ip and x denote non-relativistic quark and antiquark fields, respectively. 
The leading order Coulomb potential is part of the leading order Lagrangian 
and cannot be treated as a perturbative interaction term. The effective La- 
grangian yields the Schrodinger equation for the ti system, the only pecu- 
liarity being that the Schrodinger equation contains a non-hermitian term to 
account for top quark decay and a residual mass term Smt(nf) defined by 
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Figure 10: The total normalized ti cross section (virtual photon contribution only) 
in the threshold region at leading order (dotted), next-to- leading order (dashed) 
and next-to-next-to-leading order (NNLO) versus the center-of-mass energy.LJ Both 
plots use m t (m t ) = 165 GeV, Tt = 1.43 GeV and a s (niz) = 0.119 as in- 
put. The three curves for each case refer to the renormalization scale fi = 
{15(upper); 30(central); 60(lower)} GeV. The PS mass corresponding to fht(rn t ) = 
165 GeV is m t , PS (20GeV) = 173.30 GeV, the pole mass m t , po i c = 175.05 GeV. The 
upper panel shows the successive approximations in the PS mass scheme, the lower 
panel in the pole mass scheme. 
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Eq. (73). The residual mass implies that the non-relativistic static energy is 
\J~cj*' — 2mt l ps(p/) rather than \fq 1 ~ — 2mt and by expressing the cross section in 



terms of \/q 2 — 2m tj ps(/i/), the convergence of the perturbative approximation 
is also improved here. 

TVip l rrqss action near threshold is now known to next-to-next-to-leading 
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and it is at this order that the improv 
of the peak of the line shape becomes particularly visiblej£' 
Fig. [lfj. Subtracted top quark masses are therefore essential to take profit from 
the small experimental error. As shown above the subtracted masses can be 
accurately converted to the.MS mass which should then be a useful standard 
Subtractions different 
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from the potential subtraction can be conceived 
and are useful as far as they eliminate the leading infrared sensitive-tepa. from 
the pole mass. They have also been used in top quark production.t2irt2£ 



3.4 Synopsis of Other Results 



The analysis of renormalon divergences can be used to study non-perturbative 
corrections to many other hard QCD processes. The following gives an incom- 
plete discussion of some of these processes. 

Fragmentation in e + e~ annihilation. Inclusive single particle production 
in e + e~ annihilation i.s the, time-like analogue of deep inelastic scattering. A 
renormalon analysis£aii23'c3 predicts the leading power corrections to the dif- 
ferential cross section to be of order (A/Q) 2 , in agreement with the light-cone 
expansion, 
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and suggests the following parametrization: 



dx 



■(x,Q 2 




5(1 -z) + 



D q (x/z,y) 



l-z 



D g (x/z,n) 



(78) 



da 



L+T 



dx 



A 2 [ x dz 



[1 - z]+ 2 



D q (x/z,y) 



+ 



Dg(x/Z,fi) }, 



(79) 



where Di denotes the leading-twist fragmentation function for parton i to 
decay into any hadron, "L+T" the sum of longitudinal and transverse frag- 
mentation cross sections and the plus distribution is defined as usual. The 
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power corrections are added to the leading-twist cross sections as 



dap 
dx 



(x,Q 2 ) 



da\ 



dx 



<x,Q 2 ) + 



da\ 



dx 



-(x,Q 2 ). 



(80) 



The constants and d are to be fitted to data and depend on the order of 
perturbation theory and factorization scale \x adopted for the leading-twist 
prediction. 

Owing to energy conservation, the parton fragmentation functions disap- 
pear from the second moments 



H 



dx x 



dcr L,T 

dx 



(81) 



which can therefore be calculated in perturbation theory up to power cor- 
rections. The power expansion of the fragmentation cross section has strong 
soft-gluon singularities and the expansion parameter relevant at small x is 
A 2 / \Q 2 x 2 ). This can be related to the fact that in perturbation theory the 
hard scale relevant to gluon fragmentation is not Q, but the energy Qx of the 
fragmenting gluon. It was notedE£ao that these strong singularities lead to a 
linear A/Q non-perturbative correction to the ap and ap- This can be seen 
from 

A 
Q 

A/Q 





" A 2 " 


/ dx - x 


1 2 


Q 2 x\ 



(82) 



for any n, which also tells us that the correct 1/Q power correction is obtained 
only after resumming the power expansion at definite x to all orders. To the 
2-loop accuracy, the IR contribution to a^ appears to be related to that for 
the C-parameter.H 

The total cross section in e + e~ annihilation into hadrons is given by the 
sum of the transverse and longitudinal cross section. In this sum all power 
corrections of order 1/Q, 1/Q 2 and 1/Q 3 cancel as expected from the operator 
product expansion. 

Drell- Yan production. Drell-Yan production of a lepton pair or a massive 
vector boson, A + B — ► {7*, W, Z}(Q)+X , where X is any hadronic final state, 
presents the best studied case of a hard process with two disparate hard scales 
for which logarithmically enhanced contributions due to soft-gluon e missi on 
can be resummed systematically to all orders of perturbation theory JH9'EIiJ In 
Mellin space, the cross section dauy/dQ 2 factorises into a product of parton 
distribution functions and the quark-antiquark scattering cross section that 



43 



takes the form 



cj qq (N,a 8 (Q)) = H(a s (Q)) exp [E(N, a s (Q))] + R(N, a s (Q)), (83) 

where R(N,a s (Q)) vanishes as N — > oo, H(a s (Q)) is independent of N, and 
the exponent is given by 

rl Z N-1 _ 1 f r^^-^dk 2 
E(N,a s (Q)) = / dz— 2/ ^r[a s (fc t 2 )] 



+ B[a fl ((l-z)Q 2 )]+C[a a ((l-zW)] }. (84) 



A similar expression is valid for the thrust distribution, cf. Eq. (g5|), in which 
case we argued that it implies a linear in 1/Q non-perturbative correction 
because of the corresponding IR sensitivity of the first term involving the 
integral of the QCD coupling over low scales. For Drell-Yan production it hap 
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pens, however, that the infrared renormalon contributions arising in this way 
are exactly canceled by the divergent perturbative expansion of the function 
C[a s ]£!l The cancelation was shown explicitly in the large-/?o approximation 
and can be interpreted!.] as the cancelation of 0(1/ Q) corrections between soft- 
gluon emission at different angles. In particular, the large-angle, non-collinear 
emission is important. This conclusion is general and extends beyond Drell- 
Yan process. A recent reanalysialla of soft-gluon resummation in Drell-Yan 
production also arrives at the cancelation of the leadiiig-L/D IR contributions. 
The absense of 1/Q power corrections has been putjli-3't!^'llia into the more gen- 
eral context of Kinoshita-Lee-Nauenberg (KLN) cancelations. Knowing that 
any potential 1/Q correction would come from soft particles, but not collinear 
particles, the KLN transition amplitude can be constructed, which includes 
a sum over soft initial and final particles degenerate with the annihilating qq 
pair. The KLN transition amplitudes have no 1/fco (where ko stands for the 
energies of the soft particles) contributions (collinear factorization is implicitly 
assumed). As a consequence, the amplitude squared, integrated unweighted 
over all phase space, is proportional to dkoko, which by power counting implies 
at most 1/Q 2 power corrections. To make connection with a physical process, 
one has to demonstrate that the sum over degenerate initial states can actually 
be dispensed of. This can be shown in an abelian theory using Low's theorem. 
The generalization to QCD is still an open problem. 

Hard exclusive reactions. The theory of hard exclusive scattering is much 
less developed compared to inclusive processes. Also from the experimental 
side there is conflicting evidence. In this situation simple estimates using 
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renormalons can provide important insight. In a generic hard exclusive pro- 
cess one finds two sources of renormalon divergence and power corrections. The 
first is power corrections in the hard coefficient function, which are present in- 
dependently of the form of the hadron wave function. These correspond to 
higher-twist corrections in the hard scattering formalism. Additional power 
corrections are generated after integrating with the hadron wave function over 
the parton momentum fractions and these depend on the details of the wave 
function. These power corrections arise from the region of small parton mo- 
mentum fraction and can be associated with power corrections due to the "s oft" 
or "Feynman" mechanism for exclusive scattering. For the-simplest reaction 
7* + 7 — > 7T° and the deeply virtual Compton scatteringtl? 
both power corrections are of order 1/Q 2 . Another interesting application 



117 



1+rA 

119 



concerns the structure of the light-cone expansion of a non-local operator sand- 
wiched between vacuum and the pion that is used to define the pion distribution 
amplitude 4>-n{u): 



(0\d(0) W u(x)\7r + (p)) 



if*p M due- lupx 
Jo 



1 



-x 2 9l {u)+0{x A ) 



(85) 



The function g\(u) is interpreted as a two-particle pion distribution function 
of twist-4tH3 and is usually estimated to be gi(u) ~ const • u 2 (l — u) 2 by 
the contribution of the lowest conformal operator. The renormalon estimates 

gh orders in the conformal 



no 



g[ en (u) ~ const ■ u(l — u). 



can give constraints on possible contributions o ; 
expansion and in the large-/3o limit the result isli 
This result is significant since the behavior of gi(u) at the end points determines 
the parametric size of power corrections to inclusive reactions involving pions. 
In particular, a nonzero value gi(u — > 1) would invalidate factorization for the 
pion form factor. The further softening of the endpoint behavior (1 — u) — > 



(1- 



predicted by the conformal expansion must be interpreted as the effect 



of the hierarchy of the anomalous dimensions of conformal operators. 



Deep inelastic scattering in the small-x limit. Exploratory studies 
have been performed in the context of small-x structure functions to use renor- 
malons in order to obtain a certain prescription to implement the running 
coupling i n th e BFKL equation. An apparent 1/Q correction to the kernel is 
suppresse d! 122 ! after convolution with the hadron wave function such that the 
correction to the structure function is only of order 1/Q 2 . More recently, sim- 
ilar studios have been carried out involving the next-to-leading order BFKL 
kernel.l 123 " 124 ! The series expansion of the solution to the BFKL equation with 
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the exact 1-loop running coupling produces a series expansion of the form 

]T (ayal) n/2 T(n/2), (86) 

n 

where a = 42£(3)/3q /tt and y is the (large) rapidity that characterizes a scat- 
tering process in the BFKL limit. If we take the Borel transform with respect 
to ttg, the above series leads to a typical renormalon pole. The unusual feature 
is that location of the renormalon pole depends on the kinematic variable y, 
and not only in overall prefactor. When aya^ ~ 1 the series diverges from 
the outset and no perturbative approximation is possible. This leads to the 
constraint y < l/(aa^) for rapidities to which the BFKL treatment can be 
applied, in agreement with other methods. 



4 Ultraviolet Renormalons 



Perturbative expansions in renormalizable field theories exhibit another kind 
of divergence known as the ultraviolet (UV) renormalon. This divergence is 
easy to identify in simple sets of diagrams. For example, the integral in Eq. (^|) 
has a contribution of order (—0o/q) n nl from large k 3> Q, where q depends on 
the behavior of I(k) at large k. For UV finite or renormalized quantities the 
leading term I(k) ~ \/k A is absent or subtracted and I(k) ~ 1/fc 6 in general. 
This leads to q = 1. Higher orders in the large-fc expansion of I(k) give rise 
to divergent behavior with larger q and the corresponding singularities in the 
Borel transform. 

We note that since k^> Q the powers of the logarithm in Eq. (|J) alternate 
in sign. Hence the terms of the divergent series are also sign-alternating. For 
such series the Borel integral can be unambiguously defined and so it seems that 
we need not concern ourselves further with this type of divergence. However, 
it is not obvious that Borel summation is the correct way to recover the non- 
perturbative result. The series alone tells us only that we can approximate the 
exact result up to an uncertainty of order (A/Q) 2 , if q = 1 and if the series 
expansion is expressed in terms of a s (Q). 

This is a puzzling result. For the GLS sum rule /glSj used as an example in 
the introduction, this implies an ambiguity from arbitrarily large loop momenta 
as large as the leading higher-twist correction, even though the strong coupling 
is small in the ultraviolet. And how could such an UV effect be non-local (non- 
polynomial) in the external momentum Q? We shall see later in this section 
that the estimate (A/Q) 2 is too naive and that there is indeed no ambiguity 
associated with UV renormalons. But we shall first discuss the theory of UV 
renormalons independent of particular sets of diagrams. 
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UV renormalons are related to internal loop momenta larger than any ex- 
ternal momentum; they disappear if Green functions are computed with a UV 
cutoff not taken to infinity. This suggests that the methods of renormalization 
theory can be used to characterize UV renormalons. The simple example of 
the function I{k) discussed above shows that the first UV renormalon (q = 1) 
is related to dimension-6 operators, since this is what the first non- vanishing 
term in the large-fc expansion would yield. The first UV renormalon is partic- 
ularly simple since for dimensional reasons we can allow only a single insertion 
of a dimcnsion-6 operator. Parisfl has used this insight to compute the nature 
of the leading UV renormalon singularity in the scalar (/) 4 -theory; .-the con- 



nection with effective theories at finite UV ■cutofLwa.s also pursued.tHZI Later 
particular sets of diagrams were investigatedlH§P 9 !l 130 H 131 l and the dimension-6 
operators corresponding to them were identified. In particular, the dominant 



role of four-fermion operators was emphasized by Vainshtein and Zakharov 

Parisi's method can be usedLto determine the UV renormalon singularity in 
QCD up to an overall constantEj by means of renormalization group equations. 
These equations are formulated most easily for ambiguities or, equivalcntly, 
imaginary parts of Borel-type integrals. In QCD UV renormalons lie on the 
negative Borel axis and do not lead to ambiguities. It is technically convenient 
to consider instead the integral 

—tc+ie 

I[R](a s ) = [ die-*' 01 * B[R](t) < — < i c < |-, (87) 
J Po Po 

O+ie 

given a series expansion R in a s = a s (fi) and its Borel transform B[R]{t). The 
integral is complex and its imaginary part is unambiguously related to the 
first UV renormalon singularity at t = — 1/ (3q or the corresponding large-order 
behavior. 

The UV factorization properties of Green functions imply that the imag- 
inary part of I[R] owing to the leading UV renormalon can be represented 
as 

lmI[R]{a s ,p k ) = irC i (a s )i? 0j (a S) p fe ). (88) 

Li ' 

i 

In this equation Oi denote dimension-6 operators and i?o i the Green function 
from which R is derived with a single zero- momentum insertion of 0j. Ci{a s ) 
are the coefficient functions, which are independent of any external momentum 
Pk of R and in fact independent of the quantity R. They play the same role 
as the universal renormalization constants in ordinary renormalization. The 
coefficient function being universal, the dependence of the UV renormalon 
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divergence on the observable R is contained in the factors . These factors 
can be computed order by order in a s by conventional methods. We will 
not prove Eq. j8q ) here, but rely on the fact that familiarity with ordinary 
renormalization will make it appear plausible. 

The dimension-6 operators may be thought of as an additional term, 



(89) 



in the QCD Lagrangian with coefficients such that for any R the imaginary 
part of I[R] is compensated by the additional contribution to R from AC 
From the requirement that AC be independent of the renormalization scale fj, 
or from a comparison of the renormalization group equations satisfied by I[R] 
and Rd it can be derived that 



da. 



I) Si 



C j (a s ) = 0, 



(90) 



where 7(a s ) is the anomalous dimension matrix of the dimension-6 operators 
defined such that the renormalized operators satisfy 



(91) 



S U ^ + Hi ) °j = °- 



The unusual "—1" in Eq. ( pp|) originates from the factor l//i 2 in Eq. fl8q) . The 
solution to the differential equation (KM) can be written as 



a(a s ) = eVOSo-Jof 1 ^ F(a s )Ei(a s ), 



where 



F(a s ) = exp 



dx 



Pox 2 



(92) 



(93) 



has a regular series expansion in a s and incorporates the effect of terms of 
higher order than /3i in the /5-function and 



Ei(a s ) = exp 



dx^l 
2(3(x) 



(94) 



takes into account the (transpose of the) anomalous dimension matrix. Thus, 
the coefficient functions are determined up to a s -independent integration con- 
stants Ci. (The lower limit ceo in Eq- ( P^ ) is arbitrary. A change of ao 
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can be compensated by adjusting the integration constants.) Because the 
a s -dependence in Eq. (|87|) translates into n-dependence of large-order behav- 
ior, we deduce that this n-dependence is completely determined. Only overall 
normalization factors related to the integration constants do not follow from 
the renormalization group equation. However, these integration constants are 
process-independent numbers; they depend only on the Lagrangian that speci- 
fies the theory. It is in this precise sense that ultraviolet renormalon divergence 
is universal. This should be contrasted to the case of IR renormalons, which 
depend not only on the Lagrangian, but also on other IR parameters relevant 
to a specific class of processes. 

We now specify a basis of dimension-6 operators. In general, we are also 
interested in processes induced by external currents. Here we consider only 
vector and axial-vector currents and we let them be flavor singlets. Thus, in 
expressions like (ipMtp), a sum over flavor, color and spinor indices is implied, 
and M is a matrix in color and spinor space, but unity in flavor space. To 
account for the external currents, two (abelian) background fields v p and a^, 
which couple to the vector and axial- vector current, are introduced. Their fields 
strengths F pv = d^v v - d v v^ and H pv = - d u a p satisfy d li F' iU = j v and 
daH^ = f A . A basis of dimension-6 operators is then given by 

0\ = (i/hvV'XV'T'V) O2 = (i/'7m75?/')(V'7 m 75V') 

3 = {i> lp T A ^){^T A ^) 4 = (^ 5 T A ^n 5 T A ^) 

5 = — j abc G A V G U B G Pfl c (95) 
9s 

O e = \(^)d v F^ Or 
9s 

9s 

where the overall factors \j g k s have been inserted for convenience. We neglected 
gauge- variant operators and operators that vanish by the equations of motion. 
We also assume that all Nf quarks are massless. Chirality then allows us to 
omit four-fermion operators of scalar, pseudo-scalar or tensor type. Diagram- 
matically, they cannot be generated in massless QCD, because the number of 
Dirac matrices on any fermion line that connects to an external fermion in a 
four-point function is always odd. 

To determine the leading asymptotic behavior we only need the leading 
term in a s of Eq. (p8|). To this order F(a s ) = 1 and only the leading order 
anomalous dimension matrix is needed to compute the factors Ei(a s ). We will 



9s 
9s 
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not present this calculation!! in detail here, but discuss briefly the structure 
of the result for the correlation function of two vector currents, II(Q). The 
leading contribution is caused by mixing of the four-quark operators into the 
operator Og. Taking into account that H(Q)o s (ct s , Q) oc aj 2 at leading order, 
we obtain 

ImI[IL](a s ,Q) cx % e W*°.) a -^+M& q + 0(as) ) (96) 

where Ai is related to the largest eigenvalue of the anomalous dimension matrix 
of the four-quark operators. Converting the previous equation into large-order 
behavior gives, taking fi = Q, 

r n " =°° K {-p ) n n\ n 2 -^/ti+^ = K (-#,)" n! n* 1 ' 59 ' 1 ' 75 - 1 ' 97 }, (97) 

for Nf — {3,4,5}. The constant K is related to the integration constants Cj 
and it is not known how to compute it. The result generalizes to a non-flavor- 
symmetric vector current with a different constant K. 

It is interesting to note that subleading corrections to the asymptotic be- 
havior can be computed without introducing further "non-perturbative" pa- 
rameters in addition to the constants Cj already present at leading order. As a 
rule, to obtain the coefficient of the l/n k correction, one needs the /^-function 
coefficients 0o, . . . ,Pk+i, the (k+ I)-loop anomalous dimension matrix and the 
/c-loop correction to Green functions with operator insertions. The renormal- 
ization group treatment can in principle be extended to the next singularity 
in the Borel plane at t = —2/ Pa, but this is far more complicated, since sin- 
gle insertions of dimension-8 operators and double insertions of dimension-6 
operators would have to be taken into account. 

Now that we know how to compute UV renormalons we return to the ques- 



tion whether there exist power corrections associated with them. Equation (96) 
shows that the power-suppressed imaginary part is, after all, polynomial in Q. 
Keeping the renormalization scale fj, and Q diffcrcpt- we then find that the 
minimal term of the UV renormalon series scales aaH3 

Q 2 A 2 

— — x [powers of a s l (98) 

rather than (A/Q) 2 . Since /j is an arbitrary parameter we are free to make it 
large. The effects of this is to delay the onset of UV renormalon divergence, 
as the minimal term of the series occurs at order n ~ l/(/3oa; s (/i)), and to 
simultaneously render the minimal term small. (If the number of known terms 
in the series is small, this is not a practical way to deal with UV renormalons. 
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In this case one can make use of conformal mappings in .the. Bore] plane to 



render the UV renormalon singularities less important l 133 !! 13 ^'! 13 ^) Essentially 
this phenomenon occurs because the natural scale of the UV renormalon di- 
vergence is much larger than Q. By adjusting /i in this way, we decrease the 
perturbative coefficients. Since the strong coupling also decreases, the UV 
renormalon divergence becomes irrelevant. This is to be contrasted with the 
case of IR renormalons in which case fj, would have to be adjusted to low val- 
ues. As the perturbative coefficients decrease, the coupling increases and both 
effects exactly compensate each other as far as the minimal term of the series 
is concerned. Note that these purely perturbative considerations do not imply 
that there cannot be any non-perturbative power correction associated with 
short distances. But they tell us that whatever these corrections are they bear 
no relation with UV renormalons. 



5 Conclusion 

During the past few years infrared renormalons have emerged as a useful tool 
for the analysis of IR contributions to Feynman diagrams beyond the leading 
power in the perturbative scale, and, by implication, for the analysis of non- 
perturbative power corrections to generic hard processes. 

In this review we have given a brief survey of the basic ideas and con- 
sidered a few selected applications of the method, which we think are the 
most relevant. An important problem that we have not discussed is the re- 
lation of renormalons to lattice calculations in general and in particular to 
non-perturbative lattice renormalization procedures for higher-twist operators. 
Possible non-perturbative corrections coming from small distances (e.g. due to 
instanton contributions) have also been left out. 

We attempted to emphasize that the subject of IR renormalons is closely 
related to the extension of QCD factorization theorems on logarithmic IR 
singularities to leading power corrections in the hard scale. This extended 
factorization is at present well understood only at the 1-loop level (e.g. for 
event shapes). An outstanding problem that remains even at this level is to 
compute the (leading) logarithmic corrections to the power behavior. In the 
few cases for which the result is known, its derivation is based on the operator 
product expansion and the calculation of the anomalous dimensions of the 
relevant higher-twist operators. The extension to event shapes would require 
an understanding of the renormalization group equations for the corresponding 
shape functions (cf. Sect. ||) that have yet to be written down. 

An important formal problem that remains is the calculation of the overall 
constants of the large-order behavior. This task cannot be addressed within 
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pprt.nrhatirvn theory itself and requires new methods. Some recent develop- 



mpntd 136 H 137 l may be interesting in this respect, but it appears unlikely that an 



exact analytic solution can be found. 
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